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Abstract 

We measure the subjet multiplicity M in jets reconstructed with a successive 
combination type of jet algorithm (/ct)- We select jets with 55 < Et < 100 
GeV and \r]\ < 0.5. We compare similar samples of jets at ^/s = 1800 and 
630 GeV. The HERWIG Monte Carlo simulation predicts that 59% of the 
jets are gluon jets at y/s = 1800 GeV, and 33% at ^/s = 630 GeV. Using this 
information, we extract the subjet multiplicity in quark (Mg) and gluon (Mg) 
jets. We also measure the ratio R = = 1.91 ± 0.04(stat)^Q^g(sys). 
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I. INTRODUCTION 



The Tevatron proton-antiproton collider is a rich environment for studying high energy 
physics. The dominant process is jet production, described in Quantum Chromodynamics 
(QCD) by scattering of the elementary quark and gluon constituents of the incoming hadron 
beams. In leading order (LO) QCD, there are two partons in the initial and final states 
of the elementary process. A jet is associated with the energy and momentum of each 
final state parton. Experimentally, however, a jet is a cluster of energy in the calorimeter. 
Understanding jet structure is the motivation for the present analysis. QCD predicts that 
gluons radiate more than quarks. Asymptotically, the ratio of objects within gluon jets to 
quark jets is expected to be in the ratio of their color charges Ca/Cf = 9/4 [|l|. 

II. THE kr JET ALGORITHM 

We define jets in the D0 detector 0] with the kx algorithm [Q]. The jet algorithm starts 
with a list of energy preclusters, formed from calorimeter cells or from particles in a Monte 
Carlo event generator. The preclusters are separated by /S.TZ = \/Arp~+~A0^ > 0.2, where 
T] and are the pseudorapidity and azimuthal angle of the preclusters. The steps of the jet 
algorithm are: 

1. For each object i in the list, define da = -E^j, where Et is the energy transverse to 

the beam. For each pair of objects, also define dy = ■min{E^^, E^ j)^yi^, where D is 
a parameter of the jet algorithm. 

2. If the minimum of all possible da and dij is a dij, then replace objects i and j by their 
4-vector sum and go to step 1. Else, the minimum is a da so remove object i from the list 
and define it to be a jet. 

3. If any objects are left in the list, go to step 1. 

The algorithm produces a list of jets, each separated by ATZ > D. For this analysis, 
D = 0.5. 

The subjet multiplicity is a natural observable of a kx jet @J^. Subjets are defined by 
rerunning the kx algorithm starting with a list of preclusters in a jet. Pairs of objects with 
the smallest dij are merged successively until all remaining dij > t/cutExijet). The resolved 
objects are called subjets, and the number of subjets within the jet is the subjet multiplicity 
M. The analysis in this article uses a single resolution parameter ycut = 10~^. 

III. JET SELECTION 

In LO QCD, the fraction of final state jets which are gluons decreases with X ~ Ej^ I Y s, 
the momentum fraction of initial state partons within the proton. For fixed E^-, the gluon 
jet fraction decreases when ^/s is decreased from 1800 GeV to 630 GeV. We define gluon 
and quark enriched jet samples with identical cuts in events at ^/s = 1800 and 630 GeV to 
reduce experimental biases and systematic effects. Of the two highest Et jets in the event, 
we select jets with 55 < Et < 100 GeV and \ri\ < 0.5. 
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IV. QUARK AND GLUON SUBJET MULTIPLICITY 



There is a simple method to extract a measurement of quark and gluon jets on a statistical 
basis, using the tools described in the previous sections. M is the subjet multiplicity in a 
mixed sample of quark and gluon jets. It may be written as a linear combination of subjet 
multiplicity in gluon and quark jets: 



M = /M, + (1 - f)M, (1) 

The coefficients are the fractions of gluon and quark jets in the sample, / and (1 — /), 
respectively. Consider Eq. (|ID for two similar samples of jets at y/s = 1800 and 630 GeV, 
assuming Mg and Mg are independent of y/s. The solutions are 

^1800j\^630 _ ^630^^^1800 
1 J1800 _ J630 ^ 



(1 - M^g°o - (1 - /^^°°) 

^9 - J1800 _ J630 

where M^^°° and M^^° are the experimental measurements in the mixed jet samples at 
y/s = 1800 and 630 GeV, and f^^^^ and /^^° are the gluon jet fractions in the two samples. 
The method relies on knowledge of the two gluon jet fractions. 



V. RESULTS 
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FIG. 1. Raw subjet multiplicity in fully simulated Monte Carlo quark and gluon jets. For 
visibility, we shift the open symbols horizontally. 



The HERWIG 5.9 Monte Carlo event generator provides an estimate of the gluon jet 
fractions. The method is tested using the detector simulation and CTEQ4M PDF. We tag 
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every selected jet in the detector as either quark or gluon by the identity of the nearer (in 
r]X(j) space) final state parton in the QCD 2-to-2 hard scatter. Fig. shows that gluon jets in 
the detector simulation have more subjets than quark jets. The tagged subjet multiplicity 
distributions are similar at the two center of mass energies, verifying the assumptions in 

We count tagged gluon jets and find = 0.59 ± 0.02 and /^^o ^ o.33 ± 0.03, where 
the uncertainties are estimated from different gluon PDF's. The nominal gluon jet fractions 
and the Monte Carlo measurements at ^/s = 1800 and 630 GeV are used in Eqs. The 
extracted quark and gluon jet distributions in Fig. |T] agree with the tagged distributions and 
demonstrate closure of the method. 
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FIG. 2. Raw subjet multiplicity in jets from D0 data at ^/s = 1800 and 630 GeV. 



Figure ^ shows the raw subjet multiplicity in D0 data at ^/s = 1800 GeV is higher than 
at ^/s = 630 GeV. This is consistent with the prediction that there are more gluon jets at 
i/i = 1800 GeV compared to i/i = 630 GeV, and gluons radiate more than quarks. The 
combination of the distributions in Fig. ^ and the gluon jet fractions gives the raw subjet 
multiplicity distributions in quark and gluon jets, according to Eqs. (0-0). 

The quark and gluon raw subjet multiplicity distributions need separate corrections for 
various detector-dependent effects. These are derived from Monte Carlo, which describes 
the raw D0 data well. Each Monte Carlo jet in the detector simulation is matched (within 
ATZ < 0.5) to a jet reconstructed from particles without the detector simulation. We tag 
detector jets as either quark or gluon, and study the subjet multiplicity in particle jets M^**^' 
vs. that in detector jets M''^*. The correction unsmears M"^*^* to give MP*'^', in bins of iVf^^*. 
Figure ^ shows the corrected subjet multiplicity is clearly larger for gluon jets compared to 
quark jets. 

The gluon jet fractions are the largest source of systematic error. We vary the gluon 
jet fractions by the uncertainties in an anti-correlated fashion at the two values of ^/s to 
measure the effect on R. The systematic errors listed in Table | are added in quadrature to 
obtain the total uncertainty in the corrected ratio R = f^^lzl = 1-91 ^ 0.04(stat)^Q^g(sys). 



7 



D0 Preliminary 

D = 0,5, y„, = 1 

55 < Et'" < 1 00 GeV 

It?'"! < 0,5 

■ Quark Jets 
□ Gluon Jets 



Subjet Multiplicity M 

FIG. 3. Corrected subjet multiplicity in quark and gluon jets, extracted from D0 data. 



TABLE I. Systematic Errors 



Source 


SR 


Gluon Jet Fraction 


+0.18 
-0.12 


Jet Et cut 


±0.12 


Detector Simulation 


±0.08 


Unsmearing 


±0.04 



VI. CONCLUSION 

We extract the ycut = 10"'^ subjet multiplicity in quark and gluon jets from measurements 
of mixed jet samples at ^/s = 1800 and 630 GeV. On a statistical level, gluon jets have more 
subjets than quark jets. We measure the ratio of additional subjets in gluon jets to quark 
jets R^i 1.9 ± 0.2. The ratio is well described by the HERWIG parton shower Monte Carlo, 
and is only shghtly smaller than the naive QCD prediction 9/4. 
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